AQ 


OT-TO BE TAKEN FROM THIS ROOM 


ft 
Ae nh) 
Bit 
byte g 


PARAL 
AME bits 
Sy STH 


j 


a8 
iar 


ee 


ae 
ane 


Hi 
Bebe NTC 


iy 
AGT Se 
Baht 


wee 
Feet 


Pea Ske 
te 
ate 
SAR\ENG 
Nahas 
me 


Wi Nea Aide 


Ex wenis 
UNIERSUTATIS 
MBERTAEASIS 


Digitized by the Internet Archive 
in 2023 with funding from 
University of Alberta Library 


https://archive.org/details/Wang1983 


J J 
at ‘ivy : 
: Pore 
J 


/ 


¥ : - i. | 5 
o ODS ew oda ee 
; (decid ‘Pm mt ity Ado : bey 
el A Lay Po: rn AY A - | 
ode Mit 
Denes maak 
it af Pet iy ote 


: ; i 
Wivele ® Ve , i a ‘ Me an 5 


pea bei 


ee 


ol 


dit al 
4 i a] 


ae 
0 ag i" 


THE UNIVERSITY OF ALBERTA 


Intersection and Minimum Distance Problems for Planar Polygons 
by 


\ U ) Cao An Wang 


A THES LS 
SUBMITTED TO THE FACULTY OF GRADUATE STUDIES AND RESEARCH 
IN PARTIAL FULFILMENT OF THE REQUIREMENTS FOR THE DEGREE 


OF" MASTER TOE (SCTENCE 


Department of Computing Science 


EDMONTON, ALBERTA 


Pao oS 


se 


Abstract 

Most of the recent results in computational geometry rely 
on the attribute of convexity to provide greater efficiency. 
Given an arbitrary non-convex polygon, one approach is to 
decompose the polygon into convex pieces, then to handle each 
piece separately. On the other hand, the monotonic property 
of a polygon is also frequently used for further improvements. 

In this thesis, several fundamental problems in 
computational geometry whose solutions are based on these two 
properties, convexity and monotonicity, are studied. The 
problems can be listed as follows: 

1. polygon intersection reporting; 

2. polygon intersection detection; 

3. finding the minimum distance between two polygons 
in the £; and Ls metrics. 

The major results of the thesis are sumarized below. 

1. O(n) optimal algorithms to detect the intersection 
and to find the minimum distance between a convex 
n-gon (a polygon with n vertices) and a simple 
non-convex n-gon. 

2. O(log nn) “optimalvaloqorithms ito “find the minimum 
distance between two separable convex n-gons. 

3. O(n) optimal algorithms to find the minimum vertex 
distance between two separable convex n-gons. 

4. An O(mlog n) algorithm to report the intersection 
of a convex n-gon and a simple non-convex m-gon. 


5. An O(mlog nm) algorithm to find the minimum vertex 
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distance between a convex n-gon and a simple 
non-convex m-gon, which are non-intersecting. 

An O(nlog n) algorithm to find the minimum 
distance between two separable simple non-convex 


polygons. 
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Chapter 1 


Introdution 


1.1 The Nature of Computational Geometry 

Geometry 1S an important branch of mathematics, which is 
used in almost all areas of modern science and technology. 

The increasing reliance on the computer for solving various 
Sclentific and technological problems forces computer 
Scientists seriously to consider the relationship of pure 
traditional geometry (or non-algorithmic formulation) and 
algorithmic computational geometry. For example, a pure 
mathematician may never worry about the number of operations 
required to find the minimum distance between two planar 
polygons which contain thousands of vertices. But in the view 
of computer scientists, it 1s the number of operations that is 
an important factor for high speed computation. As a result, 
the task of translating the classical results into algorithmic 
terms iS an important one for algorithm designers. 

The pioneer in this field is Shamos, whose work 
constitutes the basis of computational geometry and provides a 
unified approach to solving many geometric problems such as 
intersection problems;inclusion problems, closest point 
problems, convex hull problems, etc. [Shamos 1975, Shamos and 
Hoey 1976]. Many computer scientists have given much 
attention not only to the fundamental problems of 
computational geometry itself, but also to its applications. 


Many areasS involve fundamental computational geometry 
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problems. In graphics, for instance, essential operations 
such as windowing, clipping, positioning, etc. require the 
solution of hidden line or hidden face elimination problems, 
which involve computation of the intersection of two or more 
objects [Freeman and Shapira 1975]. In pattern recognition, 
emphasis is put on finding the convex hull of a set of points, 
the nearest neighbors of a point, and the diameter of an 
object [Toussaint 1980]. In robot moving problems, it is 
necessary to find the minimum distance of two objects and the 
diameter of an object [Shwartz 1981]. 

All the above applications are more or less related to 
the following three fundamental computational geometry 
problems. 

(a) Polygon intersection reporting: 
Given two polygons P and Q, report their common 
intersections, i.e., report R(P)nR(Q), where R(P) 
and R(Q) represent the regions of polygons P and 
Qe 

(Db) Polygon intersection detection: 

Given two polygons P and Q, determine whether they 
intersedtideachaobherA wade ppershRCP) nRbO)=c: 

(Cc) Minimum distance between two polygons: 

Given two polygons, find the minimum distance’ 
between them and report a pair of points which are 
Separated by this distance. 


‘The minimum distance may be in the L, or L, metrics. Refer 
COP pp. “25 
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1.2 Complexity Analysis 


1.2.1 The Model of Computation 

The model of computation used in this thesis is a random 
access machine (RAM) with real number arithmetic ability. We 
assume that in planar objects, each of the following 
operations can be done in constant time: finding the 
intersection of two line segments, computing the intersection 
angle of two line segments, finding the shortest distance 
between a point and a line segment, transforming one 
coordinate representation of a point to another, and 


determining on which side of a line a point lies. 


1.2.2 Representation of a Polygon 

In this thesis, only simple planar polygons are studied. 
[Shamos 1975]. A simple polygon is a polygon whose edges do 
not intersect themselves. The vertices of a simple polygon 
can be represented by ordered pairs, in either Cartesian 
coordinates <x,y> or Polar coordinates <r,6>. 6(p), and r(p) 
will denote the polar angle and the radius of a vertex p with 
respect to a specified origin. A polygon P of nm vertices is 
called an n-gon, and is represented by a sequence of vertices 
Povlige ets One “Any paar OL*twe tecncecuiivetvertices: (p57 D-7), 
where the indices are modulo n, forms an edge of the polygon. 
A direction upon each edge is imposed such that the interior 
of the polygon R(P) always lies to the righthand side of the 


edge. In other words, the edges around the boundary of P are 
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listedntiniciockwiseaorder.o7inforderrtolavoid confusion, * Peis 
used to denote the vertex chain, and B(P) to denote the 
sequence of edges that forms the boundary of the polygon. 
From now on, tracing along the boundary of a polygon means 
tracing the boundary in clockwise order. For convenience 
sake, P is also used to denote a convex polygon, while Q 


denotes a non-convex simple polygon. 


1.3 Outline of the Thesis 
Efficient algorithms have been developed for some 
problems. The well known results for two simple n-gons can be 
Summarized as follows. 
(a) Intersection reporting problem: 
two non-convex n-gons O((nt+k)log n)? 
[Bentley 1979] 
two convex N-gons O(n) 
[Shamos 1975] 
(b) Intersection detection problem: 
two non-convex n-gons O(nlog n) 
[Shamos and Hoey 1976] 
two convex N-gons O(log n) 
[Chazelle 1980] 
(c) Minimum Euclidean distance problem: 
two convex n-gons O(1log?n) 


[Schwartz 1981] 


where kK is the number of intersection points 
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in Chapter 2, an Ol mlogmaptatigorithm for reporting the 
common intersection of a convex n-gon and a non-convex m-gon 
is briefly outlined (intersection reporting problem). The 
algorithm is simpler than the general algorithm [Bentley 1979] 
for two non-convex polygons. An O(nt+m) optimal algorithm for 
determining whether a convex N-gon and a non-convex m-gon 
intersect each other is also presented (intersection detection 
problem). 

In Chapter 3, several optimal algorithms for finding the 
minimum distance* between two separable convex n-gons, between 
two separable simple non-convex n-gons, and for determining 
the minimum distance between a convex N-gon and a non-convex 
m-gon, which do not intersect, are presented. Then the 
minimum vertex distance problem between two planar NM-gons is 
discussed. An optimal O(n) algorithm for finding the minimum 
vertex distance between two separable convex nm-gons anda 
simple O(mlog n) algorithm for finding the minimum vertex 
distance between a convex n-gon and a non-intersecting 
non-convex m-gon are presented. Finally, the minimum distance 
of two n-gons in the L, metric is studied. 

In Chapter 4, the discussion is concluded and some open 


problems for further work proposed. 


3in the L, (Euclidean) metric, unless specified otherwise. 
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Chapter 2 
Intersection Reporting and Detection for Planar Polygons 
In recent years planar geometric problems have attracted 
considerable attention. These problems arise in a number of 
applications including pattern recognition, hidden line 
elimination, etc. [Akl 1979, Avis 1981, Bentley 1979, Chazelle 
1980, Dobkin 1976, Fisher 1980, Gindy 1981, Kirkpartrick 1982, 
beesi9/79) sullen. i977). Preparatay|97 7, «Schachter, 1978, Shamos 
BOT 52) 197265) Toussaint. .1/980 .Snydeuel 90048 neparticular ithe 
problems of detecting and finding the intersection of two 
planar polygons have been investigated by a number of 
researchers [Chazelle and Dobkin 1980, Shamos and Hoey 1976]. 
As stated in Chapter 2 these problems have been solved for the 
case of two convex polygons and detection of an intersection 
for two non-convex polygons. In the case of two convex 
polygons, the convexity property 1s used. In the case of two 
non-convex polygons the planar line sweep method is applied. 
Another solution for non-convex polygons is by decomposition 
into a set of convex parts, solving each part as a convex 
polygon. In this section, the case of a convex polygon and a 
non-convex polygon is considered; the solution is based on the 
monotonic and convex property. Formally, the above planar 
polygon problems can be represented as follows, 
(a) Polygon intersection reporting problem: 
Given a convex N-gon P and a non-convex m-gon Q, 
report their common intersections (i.e., report 


R(P)nR(Q)). 
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(b) Polygon intersection detection problem: 
Given a convex n-gon P and a non-convex m-gon Q, 
determine whether or not an intersection occurs 


(i.e., ROP)nR(Q)=¢9 ?). 


Here, two algorithms are presented (one of which is 
optimal). In Section 2.1, an O(mlog n) algorithm is briefly 
outlined for reporting the common intersection of a convex 
n-gon and a non-convex m-gon (intersection reporting problem). 
In Section 2.2, an O(n+m) optimal algorithm for determining 
whether a convex NM-gon and a non-convex m-gon intersect is 
presented (intersection detection problem). Discussion of 


these results and conclusions are given in Section 2.3. 


2.1 The Intersection Reporting Problem 

For the problem of reporting the intersection of a convex 
N-gon and a non-convex m-gon, the planar sweeping algorithm, 
which is designed for two non-convex polygons, may be applied. 
That algorithm takes O(mlog n) time [Bentley and Ottmann 
1979]. Here, a different algorithm to solve the above problem 
is presented. Although this algorithm takes the same time 
order as the planar sweeping algorithm, it is somewhat 
Simpler. The main idea of the improved method is to follow 
the boundary B(Q) sequentially and find the intersection 
point(s) Yor. B(O) with? B(P)9l Im ordermtoreffacientily €indothe 
intersection points, rays for P are drawn from any point 


Oc«R(P). Each ray originates from o and passes through a 
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vertex of P, thereby partitioning the plane into n ordered 
sectors [Shamos 1975]. The sectors can then be ordered by 
their polar angles with respect to o and a binary search can 
be applied to the sectors to determine which sector contains 
the endpoints of each edge of B(Q). Thus it can be determined 
whether or not an endpoint is in R(P) and also the existence 
of an intersection of boundary B(P) with a given edge e,; in 
B(Q). If an intersection exists, then the edge e- of B(P) 


d 
which intersects e; can be found in O(log n) time. 


g 
Definition 1: Let S be an ordered set of objects, and w be a 
mapping from S to another ordered domain. S is monotonically 
increasing (decreasing) with respect to w iff for any x,yeS, 
x<y “implies w(x)<aly) (a(x)zaly)). Let xeS, S-(x)={y | yeS 
and yexiwand Setx)=iyri! yeS@and’ yx? 2 *S tsrunimodal with 
respect, to aRiffi for some x! cSpmoe( xX) eandeS— x") "are 
non-empty, and S,(x") is monotonically increasing (decreasing) 
and S.(x") is monotonically decreasing (increasing) with 


respect to w. 


Example: Assume the points in boundary B(P) are ordered such 
that x<y if y istraced before x along Pp. “Let *P ber a "convex 
Nn-gon and OcR(P) be the origin. P 1s monotonically decreasing 
with respect to a mapping 8, where 6(x) is the polar angle of 
x relative to o. Given a infinite line e, let ad be a mapping 


which maps xeB(P) to its shortest distance* from e. Then B(P) 


‘The distance is negative if x is on the other side of e. 
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is unimodal with respect to d.. 


It is easy to see that the intersection point(s) between 
the boundary B(P) of a convex n-gon P and a line segment y can 
be found in O(log n) time by making use of the unimodal 
property between B(P) and the extended line y and applying a 


binary search on the vertices of P. 


Algorithm 1 
The entire algorithm is described informally as follows. 
Bet teat pomp: Matt GP, 2 69S4hG oF Ginter Inge Degtwospolygons which 
can be represented by two linked lists. The boundary of Q, 
B(Q), is traced sequentially and the intersection(s), if any, 
of an edge of B(Q) with B(P) is recorded. For each edge 
era (anedrids one of the following four cases must occur. 
acta le edge GxrldinG gy eisetotally inside RiP), 
(i.e., @, is one of the edges bounding the 
intersection region), then continue with 
Oe = (ys, 1 Aeey) » 
bee oliiq,,<R(P)* and q,¢R(P), then e; intersects 
an edge of B(P), say A; = (Py Dy, ) Let K be 
that intersection point. Now create a new 
node kK between Ps and Dyas and set up a 
Special pointer in the linked list from k to 
Ten * 
Cae hiaogeeRP and G7 .cRiP), then iseb up a special 


‘This can be determined easily if the rays for P are drawn and 
it has been decided in which sector the endpoints of e; lie. 
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pointer from q, to the new node, say k’ which 
has been inserted between P; and Pry 


d. If both q.,q,,, £R(P) and e; does not intersect 


ttl 


B(P), then proceed directly to e 


ee However, 
Brive intersec tems (P)Mat-tworpoints:e'say Kk, 
k’, then those two nodes are inserted into 
B(P), as in the case of (b) and (c). Again a 


Special pointer is set up from k to k’. 


In order to report R(P)nR(Q), the boundary B(P) is traced 
until a node K with a special pointer is reached. Then follow 
this pointer and switch to B(Q). B(Q) is traced until a 
Special pointer is reached, following that special pointer and 
Switching to B(P) again. In this manner, by switching between 
B(P) and B(Q), their edges can be traced. The area surrounded 
by the vertex chain loop forms one of the intersection regions 
to be reported. The above procedure is continued until all of 
the vertices of P have been examined. In order to avoid 
reporting a region more than once, all special pointers are 
removed after each has been traced. 

As the algorithm references all m edges in B(Q), n edges 
in B(P), and determines whether each point of Q is in R(P), 
the time complexity of the algorithm is O(mlog n+). 

It 18 conjectured that the intersection reporting problem 
will take at least O(mlog nm) time and that the algorithm 
presented here is asymptotically optimal. The following 


theorem illustrates a lower bound for a similar problem. 
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Theorem 1: Given a convex n-gon and a non-convex non-simple 
m-gon, the time complexity for the intersection reporting 
problem is at least O(mlog n). 

Proof: Given a line and a point, assume the only operation is 
deciding on which side of the line the point lies. It is 
well known that deciding whether a point lies inside or 
outside of a convex n-gon P is at least O(log n). It can 
be shown by information theory that the same problem with 
m points is at least of O(mlog n). The above problem, 
called the m-points inel newer pro owen is reducible to the 
intersection reporting problem in linear time by 
connecting m points together to form a non-convex 
non-simple polygon Q@. A solution to the intersection 
reporting problem can easily determine whether a point is 
in a n-gon P, thereby solving the inclusion problem. Thus 
the intersection reporting problem must be at least O(mlog 


ple f 


2.2 The Intersection Detection Problem 

The intersection detection problems for two simple 
non-convex N-gons and two convex M-gons have been studied, and 
have been solved in O(nlog nm) [Shamos and Hoey 1976] and O(log 
n) [Chazelle and Dobkin 1980] time, respectively. The 
algorithms contained in [Shamos and Hoey 1976] or in Section 
2.1 solve the intersection detection problem for a convex 
N-gon and a non-convex m-gon in O(mlog nm) time. However, the 


time complexity can be further improved since both these 
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methods include some redundancy. In order to design a more 
efficient algorithm for this problem, it iS necessary to 
preprocess B(Q), so that both B(P) and B(Q) can be scanned 


sequentially, thus solving this problem in linear time. 


Definition 2: Given a polygon Q and a point o, a particular 
point g on B(Q) is visible from o, if the line segment (o0,q) 
does not intersect another point of Q@. A visible chain V(Q,o) 
is a sequence of vertices (or edges) of @ which are all 


yisSible from oO. 


The following lemma shows that each edge of @Q need not be 
compared against every edge of P in order to determine whether 


or not the two polygons intersect. 


Lemma 1: Given two polygons P, Q, a point o, oeR(P) and a 
visible chain V(Q,o), then R(P)nR(Q)46 TfFF R(P)nV(Q,0)#6 or 
ocR(Q). 
BrooOt: sa) eel part’ is obvious. 
b) 'only if part'. Assume R(P) intersects R(Q), and 
OfR(Q). Then B(Q)nR(P)4¢. For a convex polygon, there 
must exist at least one point xeB(Q)nR(P) such that x is 
visible from o. Therefore, xeV(Q,0) implies R(P)nV(Q,0)#0 


Gea 


By the above lemma, checking whether V(Q,0)nR(P)=¢ is 


sufficient to answer the intersection detection problem. The 
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following procedure describes how to construct V(Q,o), with 


two dequestss Siand Tre (Figirdy 
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Vertices are stored in two deques S and 7, where the 
concatenation of these two deques represents the portion of 
the visible chain already formed. Initially S and /T are 
empty, and the vertices, Go,Q:,-.--,Q,, are pushed into S$ (or 
T) at F, as long as the polar angles of these vertices with 
respect to Oo are monotonically decreasing (or increasing). 
Assume q; is the first vertex that backs up, Neier 
8(q; )>(<)8(q,,,)). There are two cases to be considered. 
Case 1: The edge (q,,,q;) is totally or partially invisible 

from o (Fig. 2). The following vertices Vey Vere VK 

are then removed until vertex q,, where k>/, emerges from 
either end of the edge chain (q.,q:,...,Q,;,) and becomes 
visible from o. If q, emerges from the end associated 


‘A deque is a list of elements where insertion and deletion 
can be done at both ends. 
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WA COU Omics nOwD IN MV 2 a new Verlex Gy, tie 
UIvensectlomuOrathe: ray tromvOnto.gd,., with) (G,)4Gn) ers 
pushed into S (or 7) at F with q,. If qx emerges from the 
enduassocitated with q,, theyvertices’ gq, and gy are 
pushed into S (or T) at R. 

Pacem ne ecdgenG 72,0) 1S totally Vvistble from Oo. “vertices 
Gee Oe Oy, SEEN i weace@pusned into 7 (Or S) at F as 
long as (a) their polar angles are monotonically 


increasing (decreasing) (Fig. 3) and (b) they are visible 


from © (Fig. 4). 


In any case, vertices q.,,9,,,--.,d,, are popped from S$ 
(Omomacer stint? 1: vertex q; with j<J-7 1S5"Visipie* from Oo. 
These two subcases are considered separately. 
a. With reference to Fig. 3, a new vertex 04 , 
Puep intersection Of encaray om O cOrG. , With 


edge (D5 Digs is pushed into Stor 7 }"at' 'F. 
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be Wath reference tO Fig 4 aa. few mente tec.) , 
Che intersection ob (Osa. J aWwien (Ga rQaer )s 
PSepuSshed into .(oo SJmat_r.. 

Even though special cases are considered in the above 
discussion, they have illustrated the key idea behind the 
algorithm. The above two cases are applicable to a similar 
environment in subsequent steps of the algorithm and are 
applied repeatedly with the possibility of the roles of S and 
T being interchanged. The algorithm terminates when all 
vertices in Q@ have been examined. It is easy to see that no 
vertex or edge is examined more than twice and that at most m 


new vertices are created. Thus the following lemma holds. 


Lemma 2: The visible chain V(Q,o) can be constructed in O(m) 


time. 


In order to check whether V(Q,o) intersects P, draw rays 
tom Pp from ocR(P), partitioning the plane into nm sectors; 
determine the sector to which each vertex geV(Q,0), belongs; 
and determine whether geR(P) or V(Q,O)nR(P)=¢. By the 
definition of a visible chain and the convexity of a polygon, 
V(Q,o) and P are monotonic with respect to the mapping @. 
Thus vertices in V(Q,o) and P can be scanned sequentially, 
without backtracking, to check whether V(Q,o)nR(P)J=¢. Note 
that the order of the vertices in V(Q,o) must be reversed so 
that the vertices of V(Q,o) and P are in a clockwise order 


with respect to o. This algorithm, which is similar to the 
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algorithm for merging two ordered sets, can be described 


bellow. 


Algorithm 2: Intersection detection problem for a convex 
n-gon and a non-convex m-gon. 

Input: A non-convex polygon Q=(qorGire+++Qm,) and 
Gs CONV ex pol YEON P=CpemiPs.7 fe - Pya )- 

Output: Yes, if they intersect. No, otherwise. 


Method: 


1 choose a point OeR(P) as the origin ; 


2 if “Ock(Q) 
then return ‘yes' 
else do 
3 find and reverse the order of 
ViCQEOI= C/V eck Gee ; 
4 pe ie eis: ee OO iimke i VEOUV o VSORDS ua yt) 7 


comment the index value of p is modulo fn. 


5 1<0 ; 
6 while /<k-1 
do 
¥ Vion vee O UID=) 
then do 
Wisp e os One ehem eet hand aside, Oty tym Vey 
then return ‘'yes' 
else j<(j+1)mod n; 
end 
8 else do 
if v; is on the righthand side of Os ea pe 
then return 'yes' 
else i<it+1; 
end 
end 


10 return “no” > 
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Theorem 2: The intersection detection problem for a convex 

n-gon P and a non-convex m-gon Q can be solved in O(nt+m) time. 

Proof: Steps 2 and 3 can be done O(m) time. Step 4 takes 
O(log n) time by a binary search on 6(p;). Since V(Q,o) 
and P are monotonic with respect to the polar angle 6, 
steps 6 to 9 will be repeated no more than mtn times. 


Thus, the total time is O(mtn). OD 


Theorem 3: The intersection detection problem of a convex 
n-gon P and a non-convex n-gon Q requires at least O(n) 
operations in the worst case. 

Proof seConsiderytheacase sof i twoln-gons eP=(parDisisctedd,-)% 
Calas ,Gaet merge.) pawWith thedrovertices tperturbedaslightly 
from a unit circle and laid one after another 
alternatively as shown in Fig. 5. To be precise, @(p., )< 
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ried ina be tae andMed Nh nel 2 jay 2ereenigeThe lower 
bound can be derived from the number of comparisons 
necessary to decide whether a vertex is on the righthand 
or lefthand side of an edge. Obviously, to answer the 
intersection detection problem, every vertex in P and Q 
must be also checked, and its corresponding edge in 
Polygon @) Orr, ies,» verter sp, against, edge (gad) 
anceVentex cg ec dallct COGcmO m7 >) nam Tiere hore, she 


algorithm requires no less than (2n-3) comparisons. 0 


2.3 The Intersection Detection Problem for two Nested Polygons 
An optimal algorithm for the intersection detection 
problem of convex and non-convex polygons has been presented. 

However, if R(P)nR(Q)4¢, this algorithm will not indicate 
whether or not one polygon is totally contained in the other. 
In order to solve this problem, the previously presented 
algorithm must be modified. Instead of constructing the 
visible chain V(Q,o0), find the farthest chain FlQ,o) froma 
point o, where F(Q,o) is defined as the set of points, F(Q,o) 
=x’ exeQ) xis the-farthest point from ovat a certain 
orientation}. A chain can be represented as a set of k 
vertices. Farthest chains can be constructed in linear time 
in the same way that visible chains are constructed. In order 
to check whether R(P)cR(Q) or R(Q)cR(P), it is necessary to 
determine whether V(Q,o) and F(Q,o) are both inside or outside 


of P. Obviously, this can be done in optimal linear time. 


oa) A 
_ eee 7 

Py Ties Sicaiaphiti 
is) “0s 
sqeion be sue ads mart 
—— sp 


tt a0 ei. xo2vev. & ronged 
17 alee ae : 


— ‘v 


1th Od YlauoiWeD  .¢ 
sium. ezewe ymatde 
f 7 , eae Bexm>ed>. OR » sa 


: ‘on : 
hee seks 
Se 28n a P c neta oa eed) 4 x0 nog 

>to3e AC. ht some teniave on? seal bo 
\E-0G.) asd? eeel' on @etiups’s is 20) pie 


on rm 

HoOLcew4 eotsoeted aoitoeeiszal «* 

moat arts vod mitdivop ls Pearce 

sylog xsvndo=nom bas xevues Yo mol do 3G 

{liv misisopis eldg plattsvn ai . 1696WS 
beaisaines ¥tiesos ar noe gleg sio Jon 20 tigi ; 
Hq ylevolive’ey ars meide tg @ids svlce o2 asby0 7 
nizourtecos to heoten! ,betitbomw od sam mrs ErogL 
Ora ole teatime’ of? Gert? , Vo. DIV nieds 7 diets 
si3 te beriiebvad 16,90)4% exzedy .9 10596 

O mer? thing Jeenypas~ens ef % \Or% pod 


A 
‘7 
_ 


i 


§ es Gezneserge: sd Hes Miado A |. noizetash 
ee 
ak Sesourtenoo ed nas ented teed2i0F eapke 
er a { 
iouzwssca> ems anisds ote tecdks yew omne 
423500 ré 92 4 (924592 iO 19 iyi chan sade saa 
hs 
4 tani glad a 6 'o,014 po ° 4 : or 


ge 
2 
on | 
—_ 


Theorem 4: Given a convex n-gon and a non-convex m-gon, 


detect ing whether R(P)cR(Q) or R(Q)cR(P) takes OChna ine: 
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Chapter 3 
Minimum Distance Problem 

If there iS no intersection between a convex n-gon P and 
a non-convex m-gon Q, or between two convex (respectively, two 
non-convex simple) n-gons, then finding the minimum distance 
between them may be interesting. By modifying the algorithm 
given in Section 3.2, the minimum distance between P and Q can 
be determined in O(mlog nm) time. The minimum distance between 
a convex N-gon and a non-convex m-gon can be solved in O(n+m) 
time by a monotonic property as shown in Section 2.2. 
(detadlsaare oavensine Sectiong3thied tor Im Section 3.1.2, an 
O(log n) algorithm finding the minimum distance of two 
separable convex n-gons is presented. An O(nlog n) algorithm 
for finding the minimum distance of two separable simple 
sen ies polygons is also presented. 

If the endpoints of minimum distance must be the vertices 
of given polygons, then the problem becomes more complicated. 
in Seetion 3.2.1, an optimal algorithm finding the minimum 
vertex distance between two separable convex n-gons is 
presented. In Section 3.2.2, the minimum vertex distance 
between non-intersecting a convex M-gon and a non-convex m-gon 
is considered. Finally, in Section 3.3, these problems are 


extended to the L;, metric. 
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3.1 The Minimum Distance between the Boundaries of two 


Polygons 


Definition 3: Let d(x,y) be the Euclidean distance between two 
points x and y. Let P and Q be two polygons. The minimum 
distance between P and Q is defined as d(P,Q)=Min{d(x,y) | 
xeB(P), yeB(Q)}. The minimum distance between point x and Q 
is defined as d(x,Q)=Min{d(x,y) | yeB(Q)}. If there exist two 
points Dp, Gg, such that peP, geQ, and d(p,qiiqdGPsQ), when p, a 


are called a pair of closest points for P and Q. 


Two polygons are called separable, if there exists a straight 
line 7] with the property that every vertex of P, lies on one 


side of 7] and every point of P, lies on the other. 


3.1.1 Minimum Distance between a Convex n-gon and a Non-convex 
m-gon 
It will be shown that the minimum distance between P and 
Q is the same as the minimum distance between V(Q,o) and P, 


where OcR(P). 


Lemma 3: Let p and q be a pair of closest points for P and Q. 
If P is convex then g is always visible from any point in 
R(P). 

Proof: By contradiction. Assume q is not visible from 


a point keR(P). There must exist a point q’eQ between k 
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and q. As shown in Fig. 6, a line parallel to the line 


segment (p,q) and intersecting with the line segment is 


drawn such that it passes through q’ and intersects with 


the line segment (p,k) at p’. 


Obviously, d(q,p)>d(q’,p’). 


Since p,keR(P), by the convexity of P, p’«R(P), ata pe 


d(qg’,P), contradicting d(P,Q)=d(p,q). U 
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Fig 7 


It is easily shown that there always exists a pair of 


closest points between P and Q such that one of them is a 


vertex of Por Q. Given Lemma 3, one can conclude that 


d(P,Q)=d(P,V(Q,0) )= Min({d(p,V(Q,o0))| peP} u {dlq,P) | 
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which finding the minimum distance between V(Q,o) and P may 
take O(mlog n) time. This can occur because if one scans 
V(Q,0) sequentially, as in Section 2.2, to find d(g,P) for 
some g may take O(log nm) time using a binary search on the 
sectors of P. In order to achieve the linear time bound, one 
must establish some kind of monotonic property for the visible 
chainsV/ofaQStronm PRAT Det wees xX and y be ordered in such a 
wayethat! x>yemtr xcasptracediibe fore. y/ alongevs*oThis-if 

VECV > Vsti ee tevesuey then veryeew. teva. OVerticess of¥ Plare 
ordered similarly. One wants to construct a visible chain V 
of Q@, which is monotonically decreasing with respect to the 


mapping w whereloptmaps aS@pomnt¥xeV\ tor 1tS@closestyY point in 


p! 
BtP), (icteteritex,yveVoand’ xy, then*® ol (x/s' op lV) where a> 

Gx) 0, (VIeperI eo CaTHUSI OD (Vow=2.<0,(V.1)%1 ‘Moroever, if there 
is a vertex peB(P) such that , (vj; )<p<w,(vg,), then v- >, (p)> 


V Hence, if one can construct a visible chain V such that 


tt! ° 
it has the above monotonic property with respect to the 
mapping w,, V canfiibe\scanned/ sequentially anded(V,P) %can*be 
found in linear time. Note that the direction of V must be 
reversed. 

The extreme vertices of P, are those vertices in P, which 
intersect a bounding rectangle of any fixed orientation such 
that no more than one vertex of P is contained in each edge of 
the rectangle. Without loss the generality, if the vertices 
of P are represented as CarteSian coordinates, the set of 


vertices in P with the maximum or minimum x or y coordinates 


are taken as the extreme vertices of P. The following lemma 
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shows that V(Q,S) = n V(Q,x), the set of points in Q which are 
xeS§ 
visible from all points in S, is monotonically decreasing with 


respect to mapping w where S is the set of extreme vertices 


P i 
of PP. 


Lemma 4: V(Q,S) is monotonically decreasing with respect to w, 

Pell +S Is a Setcor ‘extreme: vertices of P. 

Proof: In order to prove that V(Q,S) is monotonically 
decreasing with respect to w,, one must show that for 
any two points 
XVeV UOMO yet X>Y then w, (x)<w, ly). LpaiSesurcucrent 
to prove that (a) edges (x,uw,(x)) and (y,wp(y)) do not 
intersect, and (b) for any x<eV(Q7S). (X, a, (X)) never 
intersects with V(Q,S). 

a. Assume (X,w,(x)) and (y,w,(y)) intersect (Fig. 
8), then either (xX, w,(y) )<d(xX,wp(x)) or d(y,wp, 
(x) )<d(y,w,(y)), contradicting that (x,, (x)) 
and (y,w,(y)) are pairs of closest points. 

b. Without loss of generality, assume C, and C, 
are two adjacent extreme points for P (Fig. 
9). Since x Us@vasible from. c, and cz, ViQ,S) 
can never lie inside triangle A,C,C.x. It 1s 
easily shown that (X, wp (Xx) ) always lies inside 


triangle A,C,C.X and therefore never 


intersects V(Q,S). 0 
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Similarly, P has the same monotonic property with respect 
tO Mapplng@eerabtoens conjectuned qthateSyconsistsyofanogmore 
than Bs epommts. peUntontunatelyesiterysenot oknowngwhether there 
exists an easy method to find them. The algorithm is 


described as follows. 


Algorithm 3: Finding the minimum distance between a 
convex n-gon and a non-convex m-gon. 
EnputecA convex n-gon P and a non-convex m-gon Q. 
Output: The minimum distance and their closest points pair. 
Method: 


1 find four extreme vertices of convex n-gon P, 
and BatuwS=teatcifcs,ce!} 


2 x<; 
ee foi sto, 49do x<V(x,C 7) gend; 
V CO PSI4X; 
commenty et /(0.S)=(V5 Va Vie 


3 find siesuchethat Pinned p\VeISD 3 moa» | 


comment the index value of p is modulo n. 
i<1; D<AVo,( Py, 1P3))3 


4 while i<k-1 do 


S if os vpsp; 
then D<Min(dlv, | esos} 
j<i+1; 
6 élse DeMintd( py, lV 7V, Dy) Dade 
J-ig*1) mod in: 
end. 


Theorem 5: Algorithm 3 gives an optimal solution for finding 
the minimum distance between a convex n-gon and a non-convex 


m-gon in O(mtn) time. 
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Proof: If consider the case shown in Theorem 1, and assume the 
basic operations are computing the minimum distance 
between a point and a line segment or between two points, 


then O(mtn) operations are necessary. 0 


3.1.2 Minimum Distance between two Separable Convex n-gons 
Consider the problem of finding the minimum distance 
between two separable convex n-gons P, and P, with ,,N, 
vertices, respectively. The best known algorithm [Schwartz 
1981] takes O(log?n) by performing double binary search on 
vertices of P, and P,. Introduced here is an algorithm 
performing a single binary search on the chains of P, and P,. 
Lemma 3 shows that if po, and p. are the pair of closest points 
of P, and P,, respectively, then p, is visible from any point 
PRC Bandeisim ilar ya. Savi Ss Toler mromeanyepoint, Of RUP,).. 
Therefore, the problem of finding the minimum distance between 
P, and P, can be reduced to the problem of finding the minimum 
distance between two visible chains V, and V, of P, and P,, 
respectively. In order to design an efficient algorithm to 
find this minimum distance between P, and P, in O(log n, + log 
n,) time, the following problems must be solved: 
nia, eandaVeemusts be found invO( log ny,tlog ner time, (Lemma 
Bis 
2. it must be shown that V, and V, have the unimodal property 
with respect to the mapping A, and A,, respectively, given 
that LA me is the minimum distance between V and x (Lemma 
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3. the decreasing or increasing directions of the minimum 
distance between V, and V, must be found in constant time 


(Lemma 7); 


In order to solve the third problem, one must determine 
in constant time for any given edge (pj +P s+) in V, whether or 
not Ap; Vi )> A piri Vis). The two extended edges, e;«V, and 
eyeV. are used to determine the increasing or decreasing 
directqon—sot V4 or Vj... If Vimiand-V> are viseibleschains from 4 
Single origin, there are situations, as shown in Fig.10, for 
which the decreasing (increasing) direction of V, (V.,) cannot 
be decided in constant time. However, if V,=V(P,,S.) and 
V.=V(P.,S,), where S, and S, are sets of extreme points for P, 
and P,, then it can be shown by Lemma 7, that the third 


problem is solvable in constant time. 
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Lemma 5: V, and V, can be found in O(log n) time. 

Proof: Consider finding V(P,,0) for ocP,. Since B(P,) is 
unimodal with respect to 6, the vertex with the maximum 
polar angle and the vertex with the minimum polar angle 
can be found in O(log n) time. Thus V(P,,0) and 


V,=V(P,,S) can be found in O(log n) time. Similarly, this 
Molds TLOL avon 
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Lemma 6: V, (V,) is unimodal with respect to a. (A,), where 
Me xed (xX, V ORR APeXD Ed (XFV E72 


Proof: This lemma is proven by showing that for xX)>x;>X,,, 
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Atl 
two cases are so Similar that only the former case is 
considered. The proof is by contradiction. Assume that 
aa Xe ice ee eon wh .CX ai heXe ee bet. VY, =o, 1%) sbesthe 
elosest, pornesol Vy of comix we Now, “Eromethe definition of 
V, and Lemma 4, y., <y; Sy; (Fig. 11), and because of the 
GOnVeX Ut yCimne elt vA Cee x jet hens thesanglesiy, 0). (= 
$0" and thegangleq<90° ji Moroever, @Aa(x )>A. (xp) implies 
thateancvemp—70emand «tb 1602  Scontragicting the facuuthat 
Pugisnconvexc 15 

Lemma 7: Given a pair or edges, (p,,p..,) and (Qs Qe ? Mal We 

and V, respectively, one can determine the result of either of 

the following two comparisons in constant time, (1) 4.(p,;) vs 

Se De (2) (Qe) eVs Ai (Q5,,). 


Proof: Let the two edges (p.,p,,,) and (9549 ) be extended and 


je t 
let kK be their intersection. Without loss of generality, 


assumes thateo, 1s —petweens > mand kK.) Now). from the 


dt! 
G6ninivevon. of V4, and V.yereterrang to Fig. “12, ‘the region 


of P, can not be below the extended line (q+ oy 


jt! 
implying that (G5 Qj ) MUStMpOIntetOetnes lett. git a 
peupenatcutan sis drawn siromep  . towthemextended line 
(qa, ? intersecting at p’, there are two cases to 
consider. 
a. q; is onmihe smugntenandssidesoren-(noteuthat py” 
can be on the righthand or lefthand side of k 


(Fig. 12)). Then because of the convexity of P,, 
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clearly Az(p;)>a7(p4;): 


b. q is on the lefthand side of p’. Let x be the 


eloseste point at Vaetrom Die * Me 


Aa +X)>d( Qs yx), then clearly 


As (Qiu Jed ( Q4, x)>d(qy so cal (q;). On the other 


hand; #25 Aan »X)<d(Q;,x), a line parallel to 


(Diu 1X) can be drawn through q; intersecting 


(XP fy eatey. WiThen@tromethesconvexity of@Pyy 


yeR(P,) and bs (Qi =A G5, Xx) >A Gs py)2A1(q;). This 


will hold even if the intersection angle a>90°. 7 


Sen 


Given Lemma 7, one 
one of the chains V, or 
repeated until only two 


by a binary search, the 


P, 


can reduce in constant time at least 
V, by half. This process can be 
edges#onemfor ache) ,-remain® @ Thus, 


minimum distance between P, and P, can 


be found in O(log n,t+log n,) time. The algorithm (called 
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Algorithm 4) is omitted. The following theorem proves the 


optimality of this algorithm. 


Theorem 6: Given two separable convex n-gons P,, P., there 
exists an asympotically optimal algorithm which can find the 
minimum distance between them in O(log nm) time. 

Proof: It remains to be shown that the minimum distance 
problem requires at least log nm comparisons. Consider a 
finitersetuof positive numbers X={xepxi;.8:5X,,pe0crf the 
values of the numbers are unimodal, (i.e., there exists an 
integer aj, Osj<n, suchithat x,p>xpuforsall 0<i<j and 
X,., <x, for all j<i<n) then a lower bound can be found 
uSing reducing this to a problem, which requires at least 
log nN comparisons by information theory. Furthermore, one 
can easily show that this problem is reducible to the 
minimum distance problem in constant time. Thus the 


minimum distance problem is at least O(log n) time. 0 


Corollary: the ]ine separating two convex n-gons can be found 


fneOClog’n) Gtime. 


3.1.3 Minimum Distance between two Nested Polygons 

Two optimal algorithms to find the minimum distance 
between two convex n-gons or between a convex N-gon anda 
non-convex m-gon have be presented. In the case of one 


polygon entirely containing the second, if one disregards the 
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fact that these two polygons do intersect each other, then the 
time complexity of finding the minimum distance between their 
boundaries is linear with respect to the number of vertices in 
the polygons. The time complexity of finding the minimum 
distance between two convex n-gons will increase from O(log 
n,tlog n,) to O(n,+n,) since the minimum distances between two 
convex n-gons are multimodel with respect to the mapping, 
which maps each vertex of one polygon to its closest point in 
the other polygon. Since they are monotonic with A Rt oe 
the polar angle with an origin inside the polygon, the 
O(n,t+n,) algorithm, finding the minimum distance by 
sequentially checking the two boundaries without backtracking, 
is not difficult to design. For finding the minimum distance 
between a convex n-gon P and a non-convex m-gon Q, where 
R(Q)cR(P), the time complexity will remain unchanged. This is 
easy to show by constructing a convex hull for Q, and uSing 
the same method for the case of one convex polygon including 
the second convex polygon. The endpoint of the minimum 


distance in @ must be a vertex of its convex hull. 


Theorem 7: Given two convex n-gons P,, P. such that 
R(P,)cR(P,) or (R(P,)<R(P,)) the minimum distance between them 


can besfound iin O0nat+ne times 


Theorem 8: Given a convex n-gon P and a non-convex m-gon Q 
such that R(Q)<R(P), the minimum distance between them can be 


found in O(ntm) time. 
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3.2 The Minimum Vertex Distance between two Polygons 
In this section, the studies are restricted.to the 
minimum vertex distance between two separable planar convex 


polygons. 


ego 40 setwe wel Oy', Dy ;ceee Dane) PANG. Pel Gin, Qa emis sO yt, 
and let d(P,,P,) be the minimum vertex distance between the 
two vertex sets of P, and P,, i.e d(P,,P,)=min{d(p,q)|peP, and 
qeP,." Two vertices peP, and geP, arewecalled Closest vertices 
Ctaraeene P21 ted p,d)=d(Py,Po0. 

The minimum vertex distance problem is to find the 
minimum vertex distance and its closest vertices between two 
sets of vertices. if P, and P, are considered as two 
arbitrary sets of vertices, rather than as sequences of 
vertices from two convex polygons, then by creating the 
Voronoi diagrams, one can find the minimum vertex distance 
between P, and P, in O(nlog n) time. In this section, an O(n) 
time optimal algorithm to solve the minimum vertex distance 


problem for two separable convex n-gons is presented. 


3.2.1 Minimum Vertex Distance between two Separable Convex 
n-gons 
Definition 5: A vertex chain <x,y> of P, is defined as an 
ordered set of points (XpPp Pig rere Pz oY) where 
(Dep Pyy rene Pz? ispannorderedisubsequencesotie =Coseh, spe, /; 
andixuandiyrarespointseonpthesedges bo gpry and (Py 1Pj-,), 


respectively. An edge chain (x,y) of P, is set of edges 
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represented by the vertex chain <x,y>, 1.e., 

LOxT pr) Pts ap7, [pene GR: y)}. A visible edge chain (x,y) of 
P, from g is defined as the part of the boundary of P, that is 
visible from q, i.e., {S|SeB(P,) and line segment (g,S) does 


not contain any other points in B(P,)}. 


Fig 14 


Fid-5 


Consider a convex polygon P, and a point geB(P,). In 
Fig. 13, one can see that the minimum distance (not the vertex 
distance) between g and P, is the same as the minimum distance 
between g and the visible edge chain (x,y) of P, from q. 
Moreover, as p traces along this edge chain (x,y), d(g,p) will 
decrease gradually until some point in (x,y) is reached, from 
there on d(q,p) begins to increase (i.e., d(q,p) is unimodal 
with respect to pe(x,y)). Because of this unimodal property 
and the monotonic property between two visible chains, the 
minimum distance between two convex polygons P, and P, can be 


found in less than O(n) time. However, this unimodal property 
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is no longer valid if the vertex distance between q and the 
convex polygon P,, d(q,P,), is to be considered. For example, 
ipiFig.614,1i0(qyp) asamul.timodalmiiepais ‘takenratceachfot the 
vertices of P, in sequence. Because of this, no algorithm can 
solve this minimum vertex distance problem between two convex 
polygons in less than O(n) time. In the following, an O(n) 
algorithm for solving this minimum vertex distance problem for 
two separable convex polygons is presented. 

First, partition the vertex set of polygon P, into vertex 
chains see, ahandeer,diandnpolygon Peeinte Peeeand Brie yNow 1f 
the minimum vertex distances d(P,,,P.,),d(P.,,P,) and 
d(P,.,P2.,) can be found in O(n) time, then the minimum vertex 
distance problem can be solved in O(n) time. 

The polygon P, is partitioned as follows: A vertex q,«P, 
is chosen such that if p’«B(P,) has the property that p’ is 
Eheaclosest pointeof B( Ri.) abuomigu, cine med Comnap i )=dlqarPiv) 
(note that p’ may not be a vertex of P,), then q,,,qQy, are on 
eneasidepofatherline (gi7pS)a" Lagiact jalljjthe vertices of P, 
are ‘on themsame side of (q,,7 0 @).) Then there is at least. one 
other pair of points say q,«P. and p"«B(P,), having the same 
property that d(q, Pp" )=d(q,,Ps), where all the vertices of P, 


are on the other side of (q,,p") GE lg seas. 
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The polygon P, is partitioned by these two points p’ and p”. 
The set of vertices lying between p’ and p” which face P2 are 
called P,,=(p’,p") and the remaining set of vertices of P; are 
called Pe=( pp 7  Note*thati*p ftand p” tmay Tot’ besin IP, or 
P,2 1f£ p’ and p” are not vertices of P;. Polygon P2 is 
partitioned fhy"g Stand gq” *simillariys ire. echnf=(Cq’ fq”) cand 
22=(q",G% astshownttin “Ff igiMors, 

The line segments representing the shortest distances 
from a set of vertices to P, never intersect, by lemma 4. 
After the shortest distance line segment from an arbitrary 
vertex in Pz, to B(P,) is found in O(log n) time, the remaining 
shortest distance line segments from vertices of Pz to B(P,) 
can be found in O(n) by examining the vertices in P, 
Sequentially and scanning the boundary of P, accordingly. 
Since testing whether or not all the vertices of Pz are on one 
Side of any shortest distance line can be done in constant 
time, p’ and p" can be found in no more than O(n) time. The 


Same is true for the points gq’ and g” in B(P;2). 
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Lemma 8: The minimum vertex distance problem between P,, and 


P scan besanound in.O0.Cn).t ime. 


Proof: 


For each edge in the edge chain represented by P,,, 
Sayn Die, PD wreawpenpendigularybisector D; ‘can be 
constructed in constant time (Fig. 16). These 
perpendicular bisectors for the edges in P,, partition 
the plane into a number of regions. Because of the 
convexity of P,,, it 1S easily shown that for any 
point x lying within the region bounded by 

D7, Op Prme A; (andi cprgpz..), Aexee ganustebe equal to 
d(x,p,). Note. that some of these points, say y may 
not be bounded within two perpendicular bisectors, in 
which case d(y, p,) must be equal to the distance 
between y and one of the end points of p,. For the 
firstyvertexiqgssofnbs, ehindingethenregion!toswhich: as 
belongs may take O(log nm) time. Finding the regions 
which contain the other vertices is done by processing 
each vertex sequentially and scanning the regions in 
order. Since there are at most nN vertices and two 
scannings of the regions (forward and backward), the 
minimum vertex distance and the closest vertices 


between P,, and P;-can betfound in’O(n) time. 0 


Similarly, the minimum distance problem between P,, and 


P, can be solved. The following two lemmas describe how the 


shortest vertex distance between P,, and P., is found in O(n) 


time. 
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Lemma 9: Let p,; and Pp, be two arbitrary vertices in P,. and 


qeP., where the angle <QP; P; is greater than 90°, then the 


line segment (Q,p;?) is the minimum vertex distance /]ine 


between q and the vertex chain <Py +P; >- 


Proof: 


Consider some vertex p in the vertex chain Py 1 Ps? as 


shown in Fig. 17. Because of the convexity of P,, p 


must be on the other side of line (Ds ,P,) with respect 


to q. There are two cases to consider: 


uno 


Line segment (qg,p) intersects line segment (Ps, 
pj); then d(q,p)>d(q,p,), Since angle <QP, Py is 
greater than 90°. 

Line segment (g,p) does not intersect line segment 
(Ds Py); then since geP.,., and p «P,., gq and —p 
must be on the opposite sides of line-(q,,p’). 
This implies that the angle <gp; Pp must be greater 
than or equal to 90° because of the convexity of 
P,, and the properties of the line segment 

[Guy Dem anius d(q,p)>d(q,p;). Therefore, a(q,p, } 
is indeed the shortest vertex distance between g 


and vertex chain <Py 1 Py > O 


Pp (case 1) 
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Lemma 10: The minimum vertex distance between P,2 and 
P22 can be found in O(n) time. 

Proof; 

This lemma 1S proven by construction. Given 

Pr2=( P's Peyee +s Pp 1 P") and P22 Qe ree 45 4Q"). 

First consider the quadrangle Gq; 9, PB: Py as shown in 
B1g.THise tOnelofl thenantersortangles. must. be >90°. 
Without loss of generality, assume that angle <q, P,P; 
is greater than or equal to 90°. From Lemma 9, it can 
be concluded that d(q.,p;) is the shortest vertex 
distance from Gece P,2- Let this be the initial value 
for d(P:2,P22). Now remove vertex ”; from P22, the 
new quadrangle Dy Vs P3 and Pas wheremg; 1S the vertex 
in P22 adjacent to qj; is then considered. Again one 
of these four interior angles must be larger than 90°. 
The minimum vertex distance from one of these vertices 
to its opposite vertex chain can now be determined. 
This new minimum distance is compared with the old 
d(P,2,P2), d(P12,P2) and updated, after which one of 
these vertices 1S removed. This process is repeated 
where the shortest vertex distance d(P,2,P2z) is 
updated at each step until only two vertices remain. 
Since the removed vertices are those vertices whose 
shortest distances to their opposite vertex chains 
have already been considered, the final answer for 

GCP SPP 237i Ss “ie minimum distance between P;, and 


P22. As far as the time complexity is concerned, the 
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interior angle which is > 90° can be identified and 
the corresponding vertex can be removed in constant 
time. Since there are at most O(n) vertices in P,, 
and P,,, the minimum vertex distance between P,, and 
P,. is found in O(n) time. 0 

The algorithm (Algorithm 5) is easy to design. It is not 


given here. Thus the following theorem holds. 


Theorem 9: The minimum vertex distance problem between two 


convex separated polygons can be solved in O(n) time. 


Ps 


i+1 


Theorem 10: The minimum vertex distance problem for two 
separable convex n-gons requires at least (n-1) comparisons. 
Proof: Consider the case of two convex n-gons, P,=(Po,...,P;_; 
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chatns*(pT,'::7Daq) and Coyrer: ,qal)) sare ™ parallel eto 
each*othér@and™alsovto-the™x-axiss e(rte.s al? thevpets 
have =thewsamesy-coordrnace whileeallitherq,@senave 
another y-value) Part of the x-axis is divided into 
2nN-3 equal segments. The x-coordinate of g, must lie 
within the first segment, qj in the (2i-1)st segment 
and g,-i in the (2n-3)rd segment. Similarly for the 
Po's, the x-coordinate of the p; must lie within the 
((2n-1)-2i)th segment. One can see that 
d(q,,P,)=d(q;,p,-,) and thus the minimum vertex 
distance between P, and P, must be 
Mi ntdcg- fpr) i=190n n= ye eClearly, wo, orderstortand 
the minimum vertex distance between P, and P,, (n-17) 
comparisons are required in the worst case. 1 
With the above lemma, the proposed O(n) algorithm for 
solving the minimum vertex distance is asymptotically optimal. 
The nearest neighbor problem for a convex polygon has 
been solved by Lee and Preparata [Lee and Preparata, 1977] in 


O(n). By using Lemma 10, a Simpler algorithm for solving the 


above problem with O(n) time complexity can be designed. 


3.2.2 Minimum Vertex Distance between a Convex fM-gon and a 
Non-convex m-gon 
Consider the following problem: 
Given a convex n-gon P and a non-convex m-gon Q which do 
not intersect. Find the minimum vertex distance between them 


Such that the closest vertices are visible. 


42 its 


elldw ote: 


$i eise-x 943 FW 


43 te ozantiyo0s+x are sae 


.3aenpee a3) 1th Fone ed 


\ 
<i 
'y gods Bae heads 5h ont 7h, 8 16) 
taum .4 Bre 4 posers od ann sath 
, [2 .{t-n, cyte ite piehnia 


Sart J2@:2 2033s nw ciate ats 


“tixeiqgnes saig An Asiw me 
xe 


> 6 fhegerzed esentaltl (nes tev mom nt 
ma” 7 


nop- x xeveos~a ' 


42 


In the above problem, the closest vertex in Q need not be 
visible from the extreme points in P; therefore, the problem 
is different from that in Section 3.1.1. The above problem is 
also different from that in section 3.2.1, where the closest 
vertices may not be visible from each other. In order to 
solve this problem based on the concept of visibility, the 
definition of a visible chain must be slightly modified. If 
an area is bounded by D,, ,Dj; and (p% ,Dg)5 (Dp sD ee ), then it 
is called region i, where p.,,P;,Pj; are three consecutive 
veptifes sintPyand DPD, f¥arewperpendicular -bisectons ©(refer 
to Fig. 16). The minimum vertex distance from an arbitrary 
vertex q; <Q to convex n-gon P is aq, ,p,) LE q, lies inside 
the region ji (possiblly, the minimum vertex distance is 


aq, ,Py,) (resp. Aq; Pin)) LE q, messont07. (D7), 


a 
respectively.) If Q, is the part of the polygon Q, which is 
in region 7, ViOr Py? is ftheevrsible charnroirdctrom tp, Uwhich 
is in region 7. Define Ve Tiara ae RU ive. “VIO-P)<isithe 
concatenation of all the Sidetaenvow bei, It can be proven 
that this chain V(Q,P) has the monotonic property with respect 
to mapping w, mapping each vertex of V(Q,P) to its closest 
vertex of P. It can be shown that if p. and qy are two 
closest vertices, then both p, and qe must be in the same 
region, and q; 1s visible from the vertex p,. One can now 
conclude that the closest vertex in Q is included in V(Q,P) 
since all visible points of Q from i are also included in 
V(Q,p,). A procedure constructing V(Q,P) in O(mlog n) time is 


Simidary to thateine Section. 302.. Finding the region for which 
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an arbitrary vertex q, lies is O(log n) if a binary search on 
the nm regions of P is used. 

Start at the vertex go of Q, and find the corresponding 
regions of two consecutive vertices, say Mh and] . Connect 
Go, Qi, draw a ray from 7 to qg;:, and decide the visibility of 
dameisegnentnqi, og: )Piromep, according to the\ angle of line 
(qo,Qi), line (q:1,qz2) and ray (p,,qi). The region on which 
the vertex Q,, an arbitraryevertextge ofsQpaliesyis 
detertiinedderThen thevvisrbil@tysotsline segment s(qs, ,+qe) 
from p, is found. Only two cases can occur, which are similar 
Borehese ,iniSectioness2, iFor each qemoPeQygitstakestOClogeny 
time to find its corresponding region. “Since no vertex or 
edge of Q is considered more than twice, only m vertices are 


created, thus the following lemma holds. 


Lemma 11: A visible chain V(Q,P) can be constructed in O(mlog 
n) time. 

The proof of monotonic property between V(Q,P) and P is 
Similar togbemma §4 in Section 3,-1.5 Let d(q, ,P), Gigi, Eo be 
two consecutive minimum vertex distances, then the line 
segments (q.+P,) and (Qi, Ps? do not intersect each other 
except that Pj = Py. This property holds for the case of 


ViQ,P)J), Whe rest of the proof is 


d(p; ,V(Q,P)), dp... 
omitted. 
Therefore, one can design a procedure to check the 


minimum vertex distance between V(Q,P) and P, in O(ntm) time. 
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Theorem 11: 7he minimum vertex distance between a convex n-gon 


and a non-convex m-gon can be found in time O(mlog n). 


The lower bound for finding the minimum vertex distance 
between a non-intersecting convex M-gon and a non-convex m-gon 


remains an unsolved problem. 


3.2.3 Minimum Vertex Distance for two Nested Polygons 

Consider two convex polygons as two qe sets, 1.e., 
the minimum vertex distance may cross the boundary of the 
polygons. In this case, if two convex polygons intersect each 
Other, it is not known whether the minimum vertex distance can 
be found in O(n) time. If a convex n-gon entirely contains 
another convex n-gon and the closest vertices are visible, 
then solving the minimum vertex distance problem can be based 


on following theorem. 
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Theoremyd2:ttersPrar Plebec twosvertexchainsisuch that 
R(P,/cR(P,),. The minimum vertex distance between P, and P, 
can be found in O(n) time. 
Proof:cAssumesd (pay Pe )=dip,.Faq,) 
anosatpe; Psy =d (ip, ,ge)% 1.¢ThereforemidtpywrgnJsdlpseqr?) and 
Gt peeO,Vadl pogo, © Buted(p.,g.) can not intersect 
GV DIsq;Piaswshownm PnP FigeveZ0. en ,dis.monetonickwirth 
respect to the mapping w, which maps each vertex of P, 
into its closest vertex in P,. Therefore, the minimum 
vertex distance between P, and P, can be found in O(n) 


time. 


3.3 The Minimum Distance between two Polygons in the L, Metric 
In the previous sections, the minimum distance between 
two separable convex polygons, or between a non-intersecting 
convex polygon and a non-convex polygon was considered in the 
Euclidean (i.e., the L,) metric. The minimum distance in 
metric S is defined as follows: 
Let the two points p and g in the Euclidean plane with 
coordinates (Xp,Yp) and (Xq,Yq), respectively, anda real 
number S, where 1<S<ee be given. 

The distance between p and g in the LS metric is defined as 
ds(p,q)=(|Xp-Xq|° +|Yp-Yq|> 8 and d,.(p,q)=Max( |Xp-Xq|,|Yp-Yq|). 
Here, the minimum distance problem of two separable 
convex polygons is discussed only in the L, metric since there 

does not appear to be any practical significance in the 


general S metric. The problems are presented as follows: 
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Given two separated convex n-gons P, and P, in the 
Euclidean plane, let C,, C. be the two closest points with 
minimum distance d(P,,P,). Also let R(P,), R(P,) denote the 
Interior Of polygon =r, and Py, respectively, and V, and V> the 
corresponding sets of vertices. The minimum distance of P, 
and PF, in’ the a metriciis 
aCe, 4PJ=minidbe vi Wt xeP elyePepand d(x.VInR(P, J6REPz)=0} . 
The minimum vertex distance of P, and P, in the L, metric is 
Bee Fe Jeminidiov) lixeY. fveVe 

It will be shown that the minimum distance of two convex 
polygons can be found in O(log n) time, and the minimum vertex 
distance in O(n) time. 

The definition of visibility will be slightly modified so 


Phateitiers consistant with thesl, metric. 


Definition 6: A point peP, is L, visible from qgeP, if there 
exists a pair of connected line segments ((qg,S),(S,p)) such 
that (g,S) is parallel to the x axis and (S,p) is parallel to 
the y axis or vice versa, and ((qg,S),(S,p))nR(P,)=6 except at 


the «point «p. 


3.3.1 Finding the Minimum Vertex Distance between Two 
Separable Convex n-gons in the L, Metric 
Given two separable convex n-gons, the Cartesian 
coordinates of the extreme vertices can be found in O(log n) 
by drawing lines parallel to the x or y axis, which pass these 


extreme vertices. The lines, say Y,,Y.,X,,X., partition the 
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boundaries of the polygons into at most 10 part chains since 


the boundary of a convex polygon is unimodal with respect to 


its Cartesian coordinates. The cases are shown in Fig. 21. 


Fig, 21 


Finding the minimum vertex distance of any two part 
chains between two polygons takes time linear with respect to 
the number of total vertices in these two part chains. Let y’ 
be the intersection point of parallel lines Y’s (X’S) and the 
boundary of polygons, possibly y’eV,(V.). Let S,=<x,y> denote 

Ee vertex sate chain, where y is outside the chain <x,y’> and 
is the closest vertex to the point y’. The different cases 


are shown in Fig. 22. 
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Lemma 12: The vertex part chains S,,S.,or S,,S;, which are 

between lines Y, and Y,, are monotonic with respect to the 

mapping w which maps each vertex of S. (S;) into the closest 

vertex of S,, or maps each vertex of S, into the closest 

vertex of S,(S,). 

Proof: 

ly.  Levece), jJ=da, So) Gandud Gi =i + jahnved C14#1GS 79 
be two consecutive minimum vertex distances, where /,it7 
eS, and j,jth «S,. Assume that they intersect in k, where 
intersection means the endpoints of two distance swap as 
shown in Fig. 22. Since d(j,j) is the minimum vertex 
distance, then (k, jth)>(k,j). This is contradicting to 
the assumption that d(it+7, jth) is the minimum vertex 
distance. «Therefore, GVitie js) and: dij) do mot 
intersect each other. 

2etathe coordinates /of/ vertices’ ineS; OWS7, 'S.-ate 
monotonically decreasing (increasing) with respect to the 
y coordinates due to the convexity of the polygon. Thus, 
checking the minimum vertex distance along the vertex part 
chains need never backtrack; i.e., all mappings can be 
found in the time linear to the total number of the 


vertices in the two chains due to the convexity of S’s. 0 


Lemma 13: Given that St,.S.95% and S;,S8.,S; are not inside the 
same line pair X,,X. or Y,,Y., the minimum distance of any 
pair between the vertex chains S.,S;,S. and S;,58.,S, can be 


found in time linear with respect to the number of total 
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vertices in the pair of chains. 

Proof: Omitted, since it is enough to find 
the minimum distance only from k to S’s, where k is the 
intersection point of X’s and Y’s. And the minimum 
distance from kK to S$’S can be found in linear time. Fig. 
22 illustrates this case. 

Lemmas 12 and 13 exhaust all possible cases with two 
boundaries. In each case, the minimum vertex distance along 
the vertex part chain is checked without backtracking. Asa 
result, an algorithm to compute the minimum vertex distances 
with time linear to the total number of vertices in two 
polygons can be designed. Finally, the following two theorems 


intow gels 


Theorem 12: The minimum vertex distance in the L, metric 


between two separable n-gons can be found in O(n) time. 


Theorem 14: Finding the minimum vertex distance in the L, 


metric between two separable n-gons takes at least O(n) time. 


3.3.2 Finding the Minimum Distance between two Separable 
Convex n-gons in the L,; Metric 
Consider the minimum distance between two boundaries of 
n-gons in the L, metric. The method for solving this problem 
is similar to that given in Section 3.2. Therefore, the 
proofs which are obvious are omitted. One can prove that if 
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(resp. C2) is always visible from the extreme points in P, 
(resp. P,). Hence, finding the minimum distance between two 
convex n-gons is reduced to finding the minimum distance 
between two visible chains. Let V,=V(P,, S2) (resp. V2) bea 
visible chain, where S2 (resp. S;,) is a set of extreme points 
in the x and y axes. The visible chain can be found in O(log 
n) by uSing a binary search on P,,P2 (i.e., determine the 
endpoints of visible chains). 

The unimodel property of minimum distance between two 
visible chains also holds due to the convexity of visible 
chains V,,V2. Therefore, a solution to the remaining problem 


is shown by following lemma. 


Fi9. 23 


Lemma 14: Given a pair of edges, (Op = (De Poy 1Og = Qs ye )) in 


V;, and V2, respectively, the result of either of the following 


alk ie? 


ant senate} 61 cat Ba 
ek an born " 


- 


7 ‘ 
ix@ 20 dé s er | 
; A 7 
? 1 
ia } e 
; _ 
‘ ~ * oa . ‘ 


4 6 Sa coe 
oe 


N F oe |= 
Vetatin sidle 

ata iw 

J soneteib moninin Yo yeseqevg tobemia edt 


fy ov eo @6lon cals ertgds efet 


TSN q 
| i. 


fou los © «Sabtovedt ae ania 


wy 
. site. nivel let yt Tt he : 


a os 
a. 


- 


- 


- 


on 


two comparisons can be determined in constant time, (1) A2(p,) 
VS oha D7; ) C2) 4, (q;) vs bi (Qa) + where e,«P,, ers 
Proof: The proof method is the same as lemma 7. 

By the symmetry of the two edges. Here only 6; is 

considered. Let Hleé;) represent the half space formed by 

the infinite extended edge ej. 

Noe Tea da hu ap )ney #6 as shown in Fig. 23, then by the 
convexity of P, and P,, in the case of e; above the 
boy aby eWeter tl 15 Perpendicutarstover OLlginacingmart 
Die COeiee DL) > Mert Dey) jolt) LiemcaserOL ej below 
ENewnay hi, where MPrs=perpendicidar™tore, oniginating 
Geom chen= Ary pla) M> VAn (py) bye the convexity Ofer. 


(note that e; can not lie between rays t and A because 


of visibility) ye | : 
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rays parallel to either the x or the y axis. These 
will be called extended rays. The extended rays and 
the extended edges intersect each other. Let 

(k,cle; Jamaxi(k, j(p,),(K,I( Py, ))} and 
(ble =min tke iGo, bth) GO wet ne and Lp. TGP 7.) 
be the intersection points of extended edge e, and 
extended nay. Such that el 0 Jig) GOs anesthe, closer 
pair tok, similarly for c(e,) and b(e,), and let k be 
the intersection point of the two extended edges, B be 
the intersection angle, as shown in Fig. 24. If 

Dien) <6 {Or e,c(k, ble, ), then by the convexity of P, 


iGo Vi,)>d(p; V2). ee ersicle; J ble) ), then due to 


fue 


the convexity of P, and P,, aq, Vs )>d(q Vs 


ae 
(refer to Fig. 24) If e, is outside (k,b(e,)) and e, 
is outside (k,b(e,)), then either DG Vi pdt, Vs) 
or dip, V2)>d(p; V2). Note that either €;, Or e, 
must be one of the above cases, and at least one of 


these two comparisons can be decided. 


In any case, one can determine the shorter distance. 
Therefore, an algorithm to find the minimum distance in O(log 
nm) time can be designed. Finally, the following theorem 


holds. 


Theorem 15: Given two separable convex n-gons P, P., there 


exists an optimal algorithm which can find the minimum 
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distance between them in O(log n) time. 


3.3.3 Finding the Minimum Distance between a Convex n-gon and 
a Non-convex m-gon in the L, Metric 
In order to achieve the optimal linear time bound, the 
boundary of Q is preprocessed such that the preprocessed B(Q’) 
and B(P) have the monotonic property. The preprocessing for Q 
is similar to that in Section 4.1, and the steps are described 
as follows: 

1. Find the extreme vertices of P in x,y coordinates, say 
(XieW ), UXa Voy CMe Vi Ory 2) eee etne Text reme 
vertices, draw lines paralleling either x or y axis. 
These lines partition the plane into at most 12 
regions, for each region, assign a pair of deques, as 


shown in Fig. 25. 


bide 2) 


2. Start from go, Gi, find the regions they locate; 
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connect these two vertices, and then find the 
intersection points with the four lines mentioned 
above. If these two vertices lie in different 
regions, push the line segments into their 
corresponding deques. 

jin POL anmarbitrary vertex G7, tind, the: regione at 
locates feconunect 4(qinyg, /7#and findpthe ;points 
intersecting the four lines. For each line segment 
decide its visibility from its corresponding chain of 
P. Repeat step 3 until all vertices of Q are 


examined. In order to decide the visibility of a line 


segment, consider the following three sSubcases. 
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Refer to Fig. 26a. If q becomes visible, then 
Pooneg, Of Oo co eo OeQues andr delecemrne sper t 
of the chain which previously belongs to the 
visible chain but now becomes invisible due to the 
new line segment Ge, 1g) Ccherwise, celeren mac G 
8) ye 

Refer to Fig. 26b. If q becomes invisible, i.e., 
B<0°, then delete it. Otherwise, push it into the 
deques which are assigned to this region. If the 
two endpoints in the deques belong to the lines H, 
and L,, respectively, then block the deques 
assigned to the region 5, because no vertices in 
region 5 will be closer than that in region 11 to 
Ene wenain OL ps 

Refer to Fig. 26c. If gq becomes visible from k, 
then push it into the deques, and delete the 


invisible part. 


Each vertex of B(Q) has at most 4 new vertices 


created, 


and each vertex is checked no more than twice. 


Hence the preprocessing takes only O(m) time. Because the 


preprocessed B(Q’) and B(P) are monotonic with respect to 


either the x or the y coordinates, therefore they are also 


monotonic with respect to the mapping w, which maps each 


vertex Ot BIO”) “into its closest pointe in BCP), and vice 


VET oa. 


The following two lemmas hold. 


Lemma 15: The preprocessed B(Q’) can be constructed in 
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O(m) time, B(Q’) and B(P) are monotonic with respect to 
the mapping w, which maps each vertex of B(Q’) into its 


closest point in B(P). 


Lemma 16: The minimum distance between a convex n-gon and 
a non-convex m-gon in the L, metric can be found in O(n) 


time. 


The vertex distance between them can also be found in 
O(n) time because the minimum vertex distances between 
them have the same property. Finally, the following 


theorem holds. 


Theorem 16: The minimum (vertex) distance between a convex 
n-gon and a non-convex n-gon in the L, metric can be found 
in O(n) time. 

If a convex N-gon entirely contains another convex 
n-gon, the minimum vertex distance in the L, metric can be 
found in O(n) time. The method is as follows: find the 
extreme vertices; draw the lines parallel to the x and y 
axeS; partition two polygons into at most 12 part chains. 
Since these part chains have the monotonic property with 
respect to the x and y axes, it 1S easy to design an 
algorithm with time linear with respect to the number of 
vertices in the two polygons. 

For the minimum vertex distance between a convex 


polygon and an included non-convex polygon in the metric 
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L,, such that the closest vertices must be visible each 
other, whether or not a linear time algorithm can be 


designed is unsolved. 


3.3.4 Minimum Distance between two Separable Non-convex Simple 
n-gons 

There currently exists no sub-quadratic algorithm for 
solving the minimum distance problem between two separable 
ereaen aes Simple polygons. In this section, an O(nlog m) 
algorithm to solve this problem is presented, where m is the 
number of cusp vertices [Preparata 1977] in the polygons. The 
algorithm is based on the following facts: the minimum 
distances from the vertices in one non-convex polygon to 
another separable non-convex polygon are multimodal; The 
vertex must be a cusp vertex and has 90° degree visible angle; 
If the boundary of two non-convex polygons has been 
preprocessed, then these two preprocessed boundaries have 
pseudo-monotonic property with respect to a mapping w', which 
maps each vertex of one boundary to its closest point ina 
certain part of another boundary (Pseudo-monotonic means that 
the mapping w' maps a vertex in preprocessed boundary of Q, 
(or Q.) to the local closest point in preprocessed boundary of 
OO} (ouee.) fosteac of.mappingpaaevertes of @, (cr Q;) “to the 
global closest point in preprocessed boundary of Q. (or Q,) 
without loss of the minimum distance between two Q's, where 
local closest point is a closest point from a vertex to a 


certain part of boundary bounded by two consecutive closest 
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points or its 90° degree visible angle while global closest 
point is a closest point from a vertex to whole boundary); 
Lemma 18 will show that if each vertex 1S mapped to its local 
closest point, the minimum distance between two polygons can 
not be missed; Since the endpoints of the minimum distance 
between the two polygons must be on the preprocessed 
boundaries, then check the preprocessed boundaries is 
adequate; By using this pseudo-monotonic property, one can 
apply a divide-conquer method on the two preprocessed 
boundaries to find the minimum distance between them. 

In the following, a preprocessing procedure is first 
described, then two lemmas are proved, finally the algorithm 


is briefly discussed. 


The preprocessing procedure is as follows; 

1. Find the convex hull of Q,,Q:, and draw a separating line 
L between the two convex hulls. 

2. For each vetex of Q, (respectively, Q2) draw a line 
perpendicular to L and intersecting L at a point. If the 
perpendicular line does not intersect Q, at another point, 
the vertex of Q, is called perpendicular visible from L. 
The perpendicular visible chain of Q, (Q2) from L can be 
easily decided (Refer to Fig. 27a). 

3. Construct a vertex chain based on the perpendicular 
visible chain such that each vertex in the new chain is 
cusp vertex[Preparata 1977] (Refer to Fig. 27b). 


4, Construct a new cusp vertex chain after step 2 such that 
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each cusp vertex has 90° degree visible angle by scanning 


twice of the cusp chain. 


Else a7 


Lemma 17: The minimum distance between two separable simple 
non-convex polygons is equal to the minimum distance between 
two preprocessed boundaries which can be constructed in O(n) 
time. 

Proof: (refer to Fig. 28) Assume / is not perpendicular 
invisible from lL, and d(Q, ,QW=(i,k). *Clearly,, (j,k) is less 
than (j,k) since B>90°, contradicting the assumption. Both 
step 1, 2 and 3 take O(n) time. Step 4 takes O(m) since the 
cusp chains are monotonic with respect to a mapping, which 
maps each vertex in the boundaries to its perpendicular point 
in L and the procedure equals to trianglate a monotone simple 
polygon [Garey M. etc. 1978]. Therefore, the preprocessing 


takes only O(n).0 
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Lemma 18: The preprocessed boundary of Q, (Q@.) is monotonic 
with respect to a mapping w', which maps each vertex in the 
preprocessed boundary of @, to its local closest point in the 
preprocessed boundary of Q.. 

Proof: Case 1. Refer to the Fig. 29. If f is bounded by two 
consecutive colosest pairs (k,m) and (k’,m’), then (k,j) or 
(i,im)-is~shorter thane (47°79. 

Case 2. Refer to the Fig. 30. If / is bounded by two 
neighbor vertices which form 90° degree visible angle of ji, , 
then (m’,k’) is shorter than (m’,k) or (j,k) is shorter than 


(m,k).0 
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Chapter 4 
Conclusion 

By investigating three fundamental polygon problems, this 
thesis provides some evidence to show that efficient 
algorithms for three problems can be derived if one takes 
advantage of the monotonic and convex properties of polygons. 
The monotonic Lp paeaeny plays an important role in some 
non-convex object problems. As a result, this may be a simple 
way to solve some non-convex polygon Beobrens without using 
VORONOI diagrams. Monotonicity can be used to find the 
diameter of a convex polygon (hence, a non-convex polygon by 
constructing its convex hull) or to find the nearest neighbors 
of a convex polygon, both in O(n) time, which are the same 
time complexity as those of Lee and Praparata, [1977] and 
Shamos, [1975]. 

Extending the intersection problem to higher dimensions 
is very important to linear programming. It is still an open 
problem. 

Lower bounds have great theoretical significance. The 
oracle method [Knuth 1973] is applied on the proofs of several 
lower bounds shown in this thesis. But other non-trivial 
lower bounds remain unsolved. They may require new 
techniques. 

In the minimum distance problems, the visibility property 
plays an essential role in the design of algorithms. This is 
not unexpected in view of the close relationship beween 


classical optical physics phenomena and geometric principles. 
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The minimum distance in the L, metric has practical 


meaning if one considers phenomena in an optical anisotropic 


medium. 


All results are summarized in the table. 


Objects 


Problems 


Poiygon 
Intersection 
Reporting 
Problem 


Polygon 
Intersection 
Detecting 
Problem 


; { 
Euclidean 


(Li) 


Metric 


Minimum 


Distance 
(L,) 


Problem 


Metric 


Note: * Means the algorithm is not 


OC (n+k) log n)* 
(Zentley & 
Ottmann) 


O(nliog m)* 


to be optimal. 


two Simple 
Non-Convex 
N-gons 


A Convex N-gon 
& A Non-Convex 
M-gon 
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(boundary) 
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(vertex) 
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(boundary) 
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